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Abstract. We reanalyze the strong lens modeling of the cluster of galaxies MS2137. 3-2353 using a new UBVRIJK 
data set obtained with the ESO Very Large Telescope. We infer the photometric redshifts of the two main arc 
systems which are both found to be at z — 1.6 ± 0.1. After subtraction of the central cD star light in the previous 
F702/HST imaging we found only one object lying underneath. This object has the expected properties of the 
fifth image associated to the tangential arc. It lies at the right location, shows the right orientation and has the 
expected signal-to-noise ratio. 

We improve the previous lens modelings of the central dark matter distribution of the cluster, using two density 
profiles: an isothermal model with a core, and the NFW-like model with a cusp. Without the fifth image, the arc 
properties together with the shear map profile are equally well fit by the isothermal model and by a sub-class of 
generalized-NFW mass profiles having inner slope power index in the range 0.7 < a < 1.2. Adding new constrains 
on the center lens position provided by the fifth image favors isothermal profiles that better predict the fifth 
image properties. A detailed model including nearby cluster galaxy perturbations or the effect of the stellar mass 
distribution to the total mass inward does not change our conclusions but imposes the M/Lj of the cD stellar 
component is below 10 at a 99% confidence level. 

Using our new detailed strong+weak lensing model together with Chandra X-ray data and the cD stellar com- 
ponent we finally discuss intrinsic properties of the gravitational potential. Whereas X-ray and dark matter have 
a similar orientation and ellipticity at various radius, the cD stellar isophotes are twisted by 13° ± 3°. The sub- 
arc-second azimuthal shift we observe between the radial arc position and the predictions of elliptical models 
correspond to what is expected from a mass distribution twist. This shift may result from a projection effect of 
the cD and the cluster halos, thus revealing the triaxiality of the mass components. 

Key words. Dark Matter, Galaxies: clusters: individual: MS2137. 3-2353, Gravitational Lensing 



1. Introduction 

Cosmological N-body simulations of hierarchical struc- 
tures formation in a universe dominated by collision-less 
dark matter predict universal density profiles of halos that 
can be approximated by the following distribution 
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(1) 



The early simulations of ((Navarro et al l ll997h (hereafter 
NFW) found a = 1, leading to profiles with a central cusp 
a and an asymptotic r~ 3 slope, steeper than isothermal 
(hereafter IS). 
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* Based on observations obtained at the Very Large 
Telescope (VLT) at Cerro Paranal operated by European 
Southern Observatory. 



More recently, simulations with higher mass resolution 
confirmed that the density profile Eq. can fit the dark 
matter distribution of halos, although differ ent values of a 
were obtained by various a uthors (see e.g. iGhigna et all 
1200(1: iBullock et all EorvTbl . 

While the collision- less AC DM cosmology explains ob- 
servations of the universe on large scales, two issues con- 
cerning these halos are still debated. The first one is the 
apparent excess of sub-halos predicted in numerical sim- 
ulations, compared to the number of sate l lites in halos 
aroun d normal galaxies l|Klypin et al Lll99^lMoore et all 
1999). This discrepancy may be resolved if some of the 
sub-halos never f ormed stars in the past and are therefore 
dark st ructures iBullock et all l2001at IVerde fc Jimened . 
20021). iMetcalf fc Madaul l|200l|) : iKeetonl l)2001ald) or 



Dalai fc Kochanekl l|2002f argued that we may already see 



effects of such dark halos through the perturbations they 
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induce on the magnification on the gravitational pairs of 
distant QSOs. 

The second prediction is the existence of a cuspy universal 
profile w hich cannot explain the r otation curves of dwarf 
galaxies l|Salucci fe Burkertll2000|) . If these discrepancies 
are not simply due to a resolution problem of numeri- 
cal simulations, then, as it was pointed out by several 
authors, they may i llustrate a small-scale crisis for cur- 
rent CDM models l|Navarro fc SteinmetzL kOOft) . In or- 
der to solve these issues, alternatives to pure collision- 
less cold dark matter particles, have been proposed 
l|Spergel fc SteinhardtL Eool iBode et "all 1200 lj) . Also sev- 
eral physical mechanisms which could change the inner 
slope of mass profiles, like central super-massive black 
holes (jMilosavlievic et al 120021: lHaehnelt fc Kauffmannl 
l2002ft. tidal- merging processes inward massive halos 
llMaller fc Dekel [2002) or adiab atic compression of dark 
matter can be advocated (see e.g. iBlumenthal et al.l ll986: 
iKeetonl 12001 ah. 

The demonstration that halos do follow a NFW mass 
profile over a wide range of mass scale would therefore be 
a very strong argument in favor of collision-less dark mat- 
ter particles. Unfortunately, and despite important efforts, 
there is still no conclusive evidence that observations sin- 
gle out the universal NFW-likes profile and rule out other 
models. Clusters of galaxies studies are among the most 
puzzling. In general, weak lensing analysis or X-rays emis- 
sion models show that both singular isothermal sphere 
(SIS) and NFW fit equally well their dark matter pro- 
file, but there are still contradictory results which s eem to 
rule out either NFW or IS models (see for examnlelAlleE 



TflflilTvson et alllT^lMellimllTaClowe et al.. 200f 



Willick & Padmanabhan, 2000; Clowe & Schn eideii 12001 
Arabadiis et all l2002UAthreva et all l2002|) . This degen- 
eracy is explained because most observations probe the 
density profile at intermediate radial distances, where an 
IS and a NFW profiles have a similar r~ 2 behavior. 

A promising attempt to address the cusp-core debate 
is to model gravitational lenses with multiple arcs which 
are spread at different radial distances, where the SIS and 
the NFW slopes may diff er significantly. As emphasized 
bv iMiralda-Escude] |^)95), ideal configurations are clus- 
ters with a simple geometrical structure (no clumps) and 
with the measurements of the stella r velocity dispersio n 
profile of its central galaxy (See e.g. iKelsonet all [20021. 
The MS2137. 3-2353 cluster satisfies these requirements 
and turns out to be an exceptional lensing configuration 
with several lensed images, including a demagnified one 
we find out in this work at the very center of the lens- 
ing potential. In this paper, we analyze the possibility to 
break the degeneracy between IS and NFW mass profiles 
using new data set of MS2137. 3-2353 obtained at the VLT 
and the properties of this new fifth image. 

The paper is organized as follows. In Sect. [21 we re- 
view the cluster properties after a summary on previous 
modelings that claimed for very deep photometric observa- 
tions. This section also presents the new VLT observations 
and describes the optical properties of the cluster. Sect. El 



presents the strong lensing models for softened IS ellipti- 
cal halos and NFW cuspy profiles. We discuss the global 
agreement of both approaches within the CDM paradigm 
in Sect.0] We stress the importance of the detection of the 
fifth central demagnified image of the tangential arc sys- 
tem and discuss the observational prospects for the near 
future in Sect. Throughout this paper, we assume a 
fl = 0.3, Q. A = 0.7, and H = 100 h km s" 1 kpc -1 
cosmology in which case 1" = 3.24/i _1 kpc at the cluster 
redshift z = 0.313. 



2. The MS2137. 3-2353 lens configuration 

2.1. Overview 

MS2137. 3-23 53 is a rich cD clust er of galaxies located 
at zi = 0.313 (|Stocke et all Il99l|) . The central region 
(< 4') does not show any substructures and has a regu- 
lar visible appearance, as expected for a well dynamically 
relaxed gravitational system. The discovery of a double 
arc configuration , among which was the first radial arc 
llFort et allll992|) . makes MS2137.3-2353 a perfect cluster 
for modeling, without the need for complex mass distri- 
bution. 

The lens generates a tangential arc (A01-A02, see 
FigQ associated with two other counter- images A2 and 
A4 positioned around the cD galaxy. A01 and A02 are 
twin images with reverse parity. They are two merging 
"partial" images of the source element located inside the 
tangential caustic line. The lens potential is expected to 
produce a fifth demagnified image near the center, but the 
cD galaxy brightness peak hampers its direct detection. 
In Sect. 12. 5l we investigate in more details the presence of 
a candidate and the detection probability of this fiducial 
image. 

The lens also gives rise to a radial arc Al partially 
buried beneath the stellar diffuse component of the cD. 
This arc is a s sociat ed with the elongated image A5. 
lHammer et alJ (^397) argued the diffuse object A6 near 
A5 is probably another counter-arc associated with the 
diffuse light A'l which encompasses Al. The lens con- 
figuration is shown in Fig ^ The radial arc at about 5 
arcsec together with the tangential one at 16 arc-sec al- 
ready probe the potential at two different radii and pro- 
vide a unique way to determine its slope in this region. 
Furthermore, a radial arc together with its counter-image 
gives a stronger constraint than a tangential system on 
the potential ellipticity. 

2.2. Previous lens models of MS2137. 3-2353 

This ide al con figurati on has early prom pted 
iMellier et all l|l993h and iMiralda-Escudl l|l995f) to 
show that an isothermal elliptical model with a small core 
radius (r c < 30 hr 1 kpc ) remarkably well reproduces the 
gravitational images pattern. 

Thanks to t he hig h spatial resolution of HST images, 
lHammer et al.l l|l997h were able to confirm the lens config- 
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uration described bv lMellier et alJ and to better constrain 
the location and the shape of the counter image of the 
radial arc. They derived the properties of the mass distri- 
bution, assuming a /3-model 



p(r) = p (1 + (r/r c f) 



-0 



(2) 



with /3 = 0.87 ± 0.04, r c = 2.25 ± 0.75". This model 
confirmed that arc properties observed in lensing clus- 
ters dominated by giant elliptical galaxies can be inter- 
preted with potential well centered on their brightest clus- 
ter members. This trend is indeed robust enough to be 
generalized with a fair confidence level on similar clusters. 
Hence, only small deviations around central galaxy posi- 
tions may eventually be explored. 

For all these models the average orientation and el- 
lipticity of the potential are kept unchanged with ra- 
dius and match the st ellar light halos of the cD galaxy. 
iMiralda-Escudel 1 1995[) studied the dynamical state of the 
central stellar halo and predicted their radial velocity dis- 
persion profile. Similar studies were carried out on sev- 
eral clusters of galaxies where a tight correlation is found 
between the projected dark matter (DM ) dist ribution 
and the faint stars hal o l|Kneib et allll993l Il996|) . Later, 
iMiralda-Escudl {2002) argued that the large tangential 
deviation angle between the radial image of MS2137.3- 
2353 and its opposite counter image implies the dark mat- 
ter distribution to have a large ellipticity. It is worth notic- 
ing that self-interacting dark matter models p redict cen- 
tral halos must be circular ; so lMiralda-Escuda s argument 
may rule out these particles. 

Regarding its radial dark matter profile, despite the 
tight constraint provided by the radial arc on isothermal 
models with core, alte rnative mass profile s can naturally 
explain its properties. iBartelmannl l)l996|) demonstrated 
that the radial arc in MS2137. 3-2353 is also consistent 
with a NFW profile. It can easily produce models as good 
as isothermal spheres with core radius making the radial 
arc properties of MS2137. 3-2353 less useful than previ- 
ously expected. A primary problem was the complete ig- 
norance of the arc redshifts. Models just predicted that the 
radial and tangential arcs could be at almost the same red- 
shift, if below z w 1, or both at a large redshift. However, 
any conclusions on the inner slope of the potential are 
sensitive to these redshifts. 

Besides, in order to probe cuspy profiles one need to ex- 
plore the innermost region of the lens, where a 5th demag- 
nificd image associated to a fold arc system is expected to 
form. This task requires a careful galaxy subtraction and 
an accurate lens model which can predict whether the 
differences between the 5th image properties between a 
NFW profile and an isothermal sphere are significant and 
measurable. These goals were serious limitations to previ- 
ous modelings that could use for high resolution imagery. 
Fortunately, they are no longe r restrictions when the re- 
cent observations by Chandra l|Wise fc McNamaralEooi]) 
and by the VLT (this work) are used together with HST 
data. The new constraints provided by these new data sets 



on the geometry of the baryonic and non-baryonic matter 
components and on the lensed images properties permit 
for the first time to probe the mass profile of a cluster over 
three decades in radius, i.e. from 1 kpc up to 1 Mpc. 

2.3. New insight on the light distribution 

The HST data have been obtained from the Space 
Telescope archive. They consist in 10 WFPC2 images ob- 
tained with the F702W filter 1 . The individual frames were 
stacked using the IRAF/STSDAS package, leading to a 
final exposure time of 22,000 sec. In addition, we used 
new data sets obtained during Summer 2001 with the 
VLT/FORS instrument in optical UVI bands and with 
the VLT/ISAAC instrument in J and K 2 . The FORS 
and ISAAC data have been processed at the TERAPIX 
data center 3 . Pre-calibrations, astrometric and photomet- 
ric calibrations as well as image stacking were done using 
standard CCD image processing algorithms. We also used 
the B and R images kindly provided by S. Seitz that were 
obtained by the FORS team during the 1999 and 2000 pe- 
riods. The exposure times of these data are shorter than 
our UVI and JK data, but they are still useful for the 
photometric redshift estimates. 



Table 1. A brief summary of the VLT data. The first 
column summarizes the seeing of the final stacked images, 
the total exposure time is given in the second column. 
Also given the Magnitude Zero Points (Z.P.). The B and 
R images were obtained by the FORS team (provided by 
S. Seitz). 



Filter 


Seeing(") 


Exp. time (sec) 


Z.P (mag) 


U 


0.72 


5280 


30.856 


B 


1.2 


2400 


32.888 


V 


0.64 


6900 


33.978 


R 


0.58 


300 


32.501 


I 


0.69 


12000 


33.484 


J 


0.49 


5880 


27.643 


K 


0.50 


6480 


26.797 



The MS2137. 3-2353 optical data provide the azimuthal 
stellar light distribution and show that its geometry is el- 
liptical. Its ellipticity 4 increases with radius, starting from 
an almost circular shape at the center, and reaches quickly 
a constant value of 0.30 beyond the giant tangential arc 
location (r > 15"). The position angle is PA k, (71±4)° at 
r = 15" (see Fig [21 . Assuming a fiducial mass-to-light ra- 
tio Tj = 2 and a I-band K-correction of 0.23, we evaluate 
the rest-frame I luminosity Lj = 1.9 x 10 11 hr 1 L . 



1 Program ID: 5402; PI: Gioia 

2 Program ID: 67.A-0098(A) FORS and 67.A-0098(B) 
IS AAC; PI: Mellier 

3 http://terapix.iap.fr 

4 All ellipticities discussed here are defined as e = , 
where a and b are the major and minor axes. 
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Fig. 1. Upper left panel : overview of the lens configuration. The three arcs systems {A01,A02,A2,A4}, {A1,A5} and 
{A'f,A6}. The central cD galaxy. This F702 HST field is 56 x 56 arcsec wide (i.e. 180 x 180 h~ 1 kpc). Upper (resp. 
lower) right panel : reconstruction of arcs deduced from the single component best fit IS (resp. NFW) model fsee !3.2|) . 
In these panels are reported the observed radial arc location. The small azimuthal offset is discussed in Sect. 14.21 The 
fifth demagnified image predicted by the models near the center is detailed in Fig0J Lower left panel, detail of some 
dots used for the model fitting (see Table |A~T|) . 



The early ROS AT results of iGioia et all ljl99(t and 
lEttori fc Fabianl ill 9991) and t h e rec ent Chandra obser- 
vations of IWise fc McNamaral ifeOOll) provide additional 
clues on the cluster halo. They confirm it appears as a well 
relaxed cluster. The X-isophotes are remarkably elliptical 5 
and do not show substructures. The orientation of gas is 
almost constant PA X = 58° ± 7°, (see Fig EJ. A new m- 

5 we used the task ellipse in the IRAF/STSDAS package for 
isophotal fitting. 



teresting observational feature is the global misalignment 
between the diffuse stellar component and the hot intra- 
cluster gas. It suggests that the stellar light distribution 
docs not match exactly the DM distribution. This point 
is independently confirmed by strong lensing models and 
is discussed in Sect.lO The MS2137.3-2353 radial prop- 
erties inferred from X-rays data reveal that the brightness 
profile presents a r c ~ 7" core radius a nd an asymptotic 
slope a ~ 1.17, and an index (3 ~ 0.56. lEttori fc Fabianl 
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Fig. 2. Orientation of the isophotal major axis as a func- 
tion of radius, (squares: intra-cluster gas, diamonds: stars 
in the F702 band). The horizontal line represents the av- 
erage orientation of the DM halo from modeling beyond 
~ 8". See Sect. H 



1999|) : lAllen et alJ l)2002|) modeled the X-ray emission 

0.10 - 0.15 



rmgas 



and derived a gas mass fraction /, 
depending on the inferred cosmology. In both cases, this 
value is almost constant between 30 and 300 h^ 1 kpc. Since 
the geometry of X-ray emission follows the overall poten- 
tial and represents a small and constant mass fraction, we 
will not consider separately the gas and the dark matter 
in the lens modeling of MS2137. 3-2353 in the following. 
Instead, we will simply reduce both components to an ef- 
fective dark matter halo as the sum of the gas and true 
DM model. As a prospect, a good refinement would be 
the full introduction of this component, independantly of 
the DM halo. Moreover, a fully 3D deprojected modeling 
of both the X-ray emissivity and the strong lensing arcs 
system would certainly be the next requirement for future 
modelings. 



2.4. VLT photometry & redshifts determination 

The photometric redshifts of arcs have been mea- 
sured with the h yperz software ijBolzonella et al [2000; 
IPello et all l200lh . The redshift z p h is derived from a 
comparison between the spectral energy distribution of 
galaxies inferred from the UBVRIJK photometry and 
a set of spectral templates of galaxies which are fol- 
lowed wit h look-back time accordi ng to the evolution 
mode ls of iBruzual &: Chariot] ljl993[) (see lAthreva et all 
2002 | for details ) . The validation of hyperz is discussed in 
Bolz onella et all (l200Cf) and has been already validated us- 



ing spectroscopic redshifts on many galaxy samples. With 
the UBVRIJK set of filters, it is possible to measure all 
redshifts of our selected galaxy sample lying in the range 
0.0 < z < 3.5. The expected redshift accuracy is between 
±0.05 and ±0.2, depending on the magnitude of each arc, 
which is enough to scale the convergence of a lens model. 



For each arc, the UBVRI and J K photometry was 
done as follows. We used SExtractor llBertin fc ArnoutsL 
Il996h to estimate magnitudes in 2" apertures around a 
well defined barycenter for each part of the arcs. The V 
frame is taken as the reference since arcs are significantly 
bluer than the cD light. We also tried to take the U and 
J ones to check the robustness of the method. As well, 
results are stable against variations of aperture. 

Table 2. Photometric redshifts. Uncertainties take into 
account the scatter in the best fits with different choices 
for photometric measurements (aperture size, reference fil- 
ter...). Note that HST data are not used for photometry. 



Arc 


A01 


A02 


A2 


A4 


Zph 


1.64 ±0.15 


1.6 ±0.1 


1.7 ±0.2 


1.58 ±0.2 


Arc 


Al 


A5 


A'l 


A6 


Zph 




1.60 ±0.1 




1.1 ±0.3 
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Fig. 3. Spectral Energy Distributions resulting from pho- 
tometric redshift analysis, a) Central cD galaxy at 
z=0.313 and b) arc A5 deduced to be at z ~ 1.6. One 
can see that contrast between arcs and cD is 15 times 
higher in U then in a redder filter like F702. 

For the radial arcs Al and A'l, photometry is strongly 
sensitive to the foreground cD diffuse stellar component. 
Furthermore, Al and A'l are overlapping, so no estima- 
tion of photometric redshifts are really stable for these 
objects. A better estimation of their redshift is provided 
by their counter-arcs which both are free from contamina- 
tion. Results for all multiple images systems are summa- 
rized in table 

Taking the best determination, we conclude that z s = 
1.6 ± 0.1 for the two sources responsible of the radial and 
tangential arc systems. The models detailed in Sect. 13.27^1 
explain the need for a different redshift of the source re- 
sponsible of A'l and A6 and is consistent with the photo- 
metric redshift z B ,A'i-A6 ~ 1-1. Hence, the critical density 
at the cluster redshift and with the adopted cosmology is: 

Ecrit = 3.73 x 10 9 /iM kpc~ 2 . 
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ISand et al, | ll2002t have recently reported a spectroscopic 
determination of the redshift of arcs which are both found 
to be at z — 1.501 in remarkably good agreement with 
our color determination. From a gravitational lensing and 
mass estimate point of view, the difference between their 
spectroscopic redshift and our photometric prediction is 
un-significant. The geometric efficiency term Dd s /D s is 
slowly varying at this redshift. We nevertheless use their 
redshift estimation and the following value for the critical 
density: 

E crit = 3.78 x 10 9 /iM Q kp C - 2 . 

2.5. Detection of the fifth central image 

Gravitational optics with a smooth potential and no 
central singularities predict strong magnifi cation s hould 
produce an odd numb er of lensed images feurkeL Il98lt 
ISchneider et all ll992). More generally, the location, the 
demagnffication or even the lack of the central image are in 
principle clues on the properties of the innermost density 
profile of lenses. 

In the case of MS2137. 3-2353, we expect the large 
arc AO to have a fifth demagnified counter-image. 
Unfortunately, any simple mass models of the lens config- 
uration predicts the fifth image of this fold configuration 
should lie within one or two arc-seconds from the cluster 
center, that is inside the central cD light distribution. Its 
detection is therefore uncertain and depends on its sur- 
face brightness, its size and its color with respect to the 
cD light properties. 

In order to check whether the fifth image associated to 
AO is technically detectable, we made several lens models 
using different mass profiles which all successfully repro- 
duce the tangential and radial arcs together with their 
corresponding counter-images. We predict its position r§ 
and magnification from the softened IS and the NFW pro- 
files of Sect. 13.21 They are respectively /iis^s) = 0.2 and 
MNFw( r 5) = 0.1. The signal-to-noise ratio per HST/F702 
pixel yields : 



pix 



F702 



1 + N cI >/N 5 + N sky /N 5 



1/2 



0.7 (3) 



where N$, N c y> and N sky are the number of photo-e _ 
from respectively the fifth image, the cD and the sky 
background close to r§. Taking into account the size of 
the image, we can express the signal-to-noise in terms of 
flux, (S/N)p 702 , as a function of the magnification /ifVs) 
(assuming that the magnification fi(r A2) does not much 
change with models). 



7 x 



(4) 



F702 



The expected signal-to-noise ratio is ss 3 in the IS case and 
ps 2 for the NFW profile. So, in principle, the fifth image 
of the MS2137. 3-2353 lensing configuration is detectable. 



Using the counter- image A2 of area A A2 and flux F A2 , 
we reconstructed the predicted fifth image satisfying : 



F5 = K r 5) X M r A2) 1 X F A2 

A 5 = At(r 5 ) x Li(r A2 )~ 1 x A A2 



and inserted it inside the cD galaxy at several positions 
close (but different) to the expected location r$. We then 
determined the significance of several extraction-detection 
techniques on the Space Telescope image. A Mexican-hat 
compensated filter turned out to provide the best cD light 
subtraction and an optimal detection of the fifth image 
twins we put inside at different positions. In all cases it 
was detected exactly at the right position, whatever its 
location inside the cD and with the expected signal-to- 
noise. 

Because we used a compensated filter which smoothes the 
signal, this later is not straightforward and we had to com- 
pare the amplitude of the flux contained in the extracted 
object to the variance of the background contained inside 
independent cells of similar size ranging along concentric 
annuli located at the radius where simulated fifth image 
twins are putted (0.6" < r < 0.9"). The averaged S/N 
found in annuli is 2.6, but it scatters between 1.3 and 3.5 
depending on the local noise properties. 

Table 3. Properties of the fifth images on real data (R.) 
and predicted from the best lens modeling (IS or NFW). 
The positions (x, y) are given in arc-sec, with respect to 
the cD centroid. Position angles (P.A.) are given in de- 
grees and a/b is the axis ratio. The errors are found from 
the changes when varying some SExtractor parameters. 
Although both the position angles and the ellipticity of 
the IS and NFW are compatible with the data, there is a 
significant difference in positions. The offset |<5cc| between 
the IS and the real position is only 0.16", whereas it is 
0.36" for NFW, which is larger than uncertainties on ob- 
servations (third column |<5cc| ~ 0.05") 



ID 


(x, y) 


|£x| 


P.A. 


a/b 


S/N 


R. 


(0.64;0.70) 


0.05 


28 ± 14 


3.1 ±1.3 


2.6 


IS 


(0.52;0.81) 


0.16 


15. 


2.2 


2.6 


NFW 


(0.33;0.52) 


0.36 


27. 


2.2 


2.1 



The application of the extraction technique on the real 
data is straightforward. The brightest residual in the fil- 
tered frame shown in the right panel of Fig 01 is detected 
at the expected location when compared to models and is 
clearly the most obvious object underneath the cD. The 
object properties are listed in Table |3J They are remark- 
ably similar to the IS and NFW fifth image predictions. Its 
coordinates are however closer to the IS fifth image than 
the NFW model. The signal-to-noise ratio of the candidate 
is ~ 2.5, in very good agreement with our expectations. In 
the frame of Fig ^ the centroid position of the candidate 
is at 

r 5 = (0.64"; 0.70"), r 5 = 0.9". (5) 



R. Gavazzi et al.: Radial mass profile of MS2137. 3-2353 



7 



Candida tc 


■eoT .model 


NFW model 


* J 

simulated. ■ « 

f W 


ATP 




k 
\ 




U. cD pen L ra idua 1 L r 


i ri 











Fig. 4. Visualization of residuals from subtraction of the cD brightness. Orientation is the same as in Fig^ Left panel: 
HST filtering. One can see the real residual as well as a simulated object derived from the IS model shown on the right 
panel. Note the comparable signal-to-noise. Also plotted the Mexican-hat filter size with both positive (r < 1.5 px) 
and negative areas (1.5 px < r < 3.5 px) Also shown the noise-free reconstructed 5th image of the fold arc expected 
from IS (middle panel) and NFW (right panel) strong lensing best fit models centered on the cD. The position is 
consistent with the IS model. Note that the filtering algorithm reduces the high frequencies in the upper frame. Thus, 
shape measurements just provide rough estimates of ellipticities, but accurate orientations. Note also that the isophote 
plotted here, is illustrative and does not correspond to the any quantitative shape measurement. 



Despite its poorly resolved shape, the candidate exhibits 
an orientation PA ~ 28 ± 14° and an axis ratio a/b = 
3.1 ± 1.4, in good agreement with the values predicted by 
both models (see Table OJ). It is worth noticing that even 
the morphology of the fifth image shows similarities with 
the reconstructed images. In particular, it shows a bright 
extension inward and a smaller faint spot outward as if 
it would be dominated by two sub-clumps which are also 
visible on predictions of Fig.Q] In the following, when us- 
ing the fifth image knowledge, we apply different weights 
on the detected features for the lens modeling. The bright- 
est part of the fifth image is statistically significant and is 
associated to the dot labeled (10) on Fig.QJand table I A. 1*1 
The mappings of the other labeled dots are not as signifi- 
cant in the fifth image. Thus, we apply a ~ 9 times smaller 
weight (i.e. 3 times larger errorbars). In other words, the 
fifth image is almost reduced to a point-like information 
without shape measurements. 



3. The dark matter distribution in MS2137. 3-2353 

In this section, the properties of the dark matter distri- 
bution of MS2137. 3-2353 are discussed in view of the most 
recent constraints we obtained from VLT data. We first 
revisit a single potential model using only strong lensing 
data but no fifth image. We then compare the projected 
mass profiles of the best NFW and IS models, extrapolated 
beyond the giant arcs positions, with the weak lensing 
analysis. Finally the fifth image is included in the strong 



lensing model which is used together with the weak lens- 
ing and the cD stellar halos to produce a comprehensive 
model of the different mass components. 



3.1. Strong lensing optimization method 

The opti mization have been carried out with the 
lensmodel 6 ( KejjtoJ, |2j)Qlt ) inver sio n software. This al- 
terna tive to the lMellier et alJ l|l993() or lKneib et all l|l993[ 
Il99(fl algorithms allows us to check the efficiency and 
the accuracy of this software for arc modeling and to 
take advantage of its association tool for multiple point- 
images. This facility was initially developed by Keeton for 
multiple-QSOs but turns out to be well suited for HST 
images of extended lensed objects. The images association 
is performed by identifying conjugated substructures like 
bulges in extended images. Because of the surface bright- 
ness conservation, brightest areas in an image map into 
the brightest of the associated ones. 



Our modeling started by identifying the brightest con- 
jugate knots in each image. More precisely, when the iden- 
tification of Mcnots distinct features in images is completed 
(with respectively iVj - =i..jv knots multiplicity) one can write 
A^knots times the lens equation relating source and image 



http : //astro .uchicago . edu/~ckeeton/gravlens/ 
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positions and the lens potential <f>: 

e 1 - v^(0!) 







Sj - 



e Nj -v<i>(e Nj ] 



(6) 



This yields the following \ 2 definition calculated in the 
image plane : 



J2 X/^'rV'^V 



with 



j = l t=l 

iVknots 



1)-JV, 



par 



(7) 



(8) 



degrees of freedom, where iVp ar is the number of free pa- 
rameters in the model. Here, S^ 1 is the error matrix 
for the position of knot j in the image i and 5xij = 

•Eobs.ij ^mod 

: ij . Analogous \ 2 minimization can be done 
in the source plane in order to speed up the convergence 
process. It is only an approximation of the previous one 
that does not directly handle observational errors in the 
image plane. However, it is much faster because it does 
not need to invert Eq. JBJ). Once the minimum location 
is roughly found, one can use the image plane Ximg ^° 
determine the best parameter set with a better accuracy. 

It is worth noting that the uncertainties in the conju- 
gate points positioning done during the association pro- 
cess dominates the astrometric errors in the position of 
each knot. Typically, the systematic uncertainty is of or- 
der 0.1 — 0.2". The VLT color similarities were also used to 
confirm the associations. The mapping between extended 
images is given by the magnification matrix a : 



<h i Ho) = s ij 



(9) 



where &y is the Kronecker symbol. Hence, other fainter 
conjugate dots become easier to identify once the local 
linear transformation between multiple images is known. 
The procedure "get constraints" - "fit a model" can be it- 
erated to use progressively more and more informations. 

In MS2137. 3-2353, we kept 13 unambiguous quintuple 
conjugated dots in the tangential arc A01. Each one is as- 
sociated to four different dots in A02, A2, A4 and the fifth 
image. We selected also 6 dots in the parts of A2 and A4 
that are only triply imaged. Likewise, Al is decomposed in 
two symmetric merging images and is also associated with 
the Eastern part of A5 (6 triple conjugated dots). Fig. ^ 
and table |A~T1 summarize the associations we selected. 

The various models are actually over-constrained. The 
6 free parameters are detailed in the following section. 
Following Eq. ©, the number of constraints is 

N = 2 x {13(5 - 2) + 6(3 - 2) + 7(3 - 1)} = 118. (10) 

The first term corresponds to the regions of the tangential 
system which are imaged five times, whereas the second 
term refers to regions imaged three times. The third term 
correspond to the radial system which is imaged three 



times. 7 Nevertheless, only 25 of these 118 constraints ap- 
pear significant to represent the first and second shape 
moments of arcs, the rest is for higher order moments and 
have less weight in the modeling. 

A galaxy at the eastern part of A02 should weakly per- 
turb its location and shape. This galaxy was introduced 
in previous models but turns out to have negligible con- 
sequences. Indeed, only upper limits on its mass (a v < 
150 km/s) arise when modeling. Its introduction appears 
marginally relevant for the study and is ignored hereafter 
although its effect is shown on Fig [2 

3.2. Strong lensing models without the fifth image 
3.2.1. Dark matter density profiles 

We model the dark matter halo with two different den- 
sity profiles. In order to focus on the main differences be- 
tween isothermal and NFW profiles, we keep the models 
as simple as possible and do not include peculiar features, 
like cluster galaxy perturbations. The center of potential 
is allowed to move slightly within 2 arc-sec around the cD 
of the cD galaxy. No prior assumptions are made about 
the ellipticity and the orientation of the dark matter halo 
relative to the light nor to the X-rays isophotes. 

The first profile is an elliptical isothermal distribution 
with core radius of the form, 



p{r) 



Pc 



i + 0/V c ) 2 

which is projected in, 
K (r, 9) 



in spherical approx. 



(11) 



e = r v /l-ecos(2(0-0 o )). (12) 



The core radius r c , scale parameter b = r c + \jB? e + rj? , 
ellipticity e and position angle 9o are free parameters. R e is 
the Einstein radius and b is related to the cluster velocity 
dispersion by 



b = ( — Y Dds 
~ V 186.5 km.s- 1 ) D s 



arc — sec. 



(13) 



The second profile is an elliptical NFW mass distribu- 
tion. The 3D profile has the form 

Pc$c 



p{r) 



(r/r s ){l + r/r s ) 2 



in spherical approx. (14) 



where r s is a scale radius, p c is the critical density of the 
universe at the redshift of the lens, and S c a concentration 
parameter related to the ratio c = r2oo / r s by 

„3 



200 



3 ln(l + c)-c/(l + c) 
The convergence n writes 



(15) 



n(r, 9) = 2k, 



1 - F(x) 



1 



with x = £/r s , (16) 



7 Eq. 1101 : the —2 in parenthesis includes the unknown 
source position and the unused fifth central image. When using 
this later, the number of constraints is 156. 
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where £ has the same meaning than before, k s 
p c 5 c r s /T, crit and 



F(x) = 



tan 1 y/x 2 — 1 , for x > 1 
L= tanh -1 VI -£ 2 , for a; < 1 (17) 



i 



, for x = 1 



3.2.2. Single halo best models 



The inversion leads to two models that fit the strong 
lensing observations equally well. They reproduce the 
multiply- imaged lens configuration of both radial and tan- 
gential arcs. The NFW best fit model leads to a Xnfw = 
2.42 per degree of freedom and xfs = 3.0 for the isothermal 
profile. 8 The final model parameters and errors bars are 
summarized in table 0] The centering of the dark matter 
halo relative to the cD galaxy is discussed in Sect. I3.4T1 

The associated counter-image of the radial arc Al 
(bright and thin structure) corresponds only to a small 
part of A5 that is triply imaged. Besides, the diffuse com- 
ponent A'l can be associated to A6 only if the correspond- 
ing source is at a lower redshift than arc A1-A5. This cor- 
roborates pho tometric redshi fts results and was previously 
mentioned bv lHammer et all Here, we find the source red- 
shift Zs(A6-A'i) to be 1.1 — 1.3. 

The velocity disp ersion derived for the IS model is con- 
sistent with results o flMellier et al.ll)l993|) . The core radius 
proposed by these authors is higher because of its different 
definition. They used a pseudo-isothermal projected grav- 
itational potential 9 ; instead, we directly model the cluster 
projected density profile. Nevertheless, to ensure the same 
Einstein radius with the same central velocity dispersion 
between their model and ours, the core radius they re- 
ported must be twice the one we found. Thus, core radii 
are consistent. 

3.3. Mass profile of MS2137. 3-2353 from weak lensing 
analysis: 

Although the error bars are large, the concentration 
parameter found for the best NFW model is about twice 
the expectations from numerical simulations and from 
the current measurements done in othe r clusters, even 
those with strong lensing features (See e.g. lHoekstra et all 
2002). Since clusters are believed to be triaxial, it may 
happen that the major axis lies along the line of sight, 
increasing atifi cially its concentra tion by projection ef- 
fects (See e.g. I.Iing fc SutoL 120021) . However, it is likely 
that our model also mix together the contributions to the 
effective concentration of the cluster and of the central 
cD potentials. This possibility can be tested by compar- 
ing the strong lensing model with the weak lensing anl- 
ysis that only probes the radial mass profile at larger 

8 These values are higher than 1 but we remind that the 
models are significantly over-constrained. 

9 leading to a convergence : /t(r) 
to be compared to Eq. 1121 



r c (1 + ^2)3/2 



with x = r jr c 



scale, where the cD contribution is negligible. On large 
scales, we used the VLT images to build a weak lens- 
ing catalog of background galaxies covering a 6.4' x 6.4' 
field of view. At the cluster redshift, this corresponds to 
a physical radius of 700/i~ 1 kpc. The detailed descrip- 
tion of the catalog analysis, namely PSF anisotropy cor- 
rections and detailed galaxy weighting scheme and se- 
lection, is beyond the sc ope of this work. Th e method 
we used can be found in l Athreva et, all l2~002|K Here, we 
only compare the result of the strong lensing mass pro- 
file and the fit of the azimuthally averaged shear on scales 
100 h _1 kpc < r < IMpc. 

The shear profile is determined by using a maximum- 
likelihood analysis, based on a x 2 minimization: 



E 

8=1 



g(r l ,z i )\ 7 



(18) 



where ej is the complex image ellipticity, g the complex 
reduced shear, the photo- z and o e ^ the dispersion com- 
ing from both the intrinsic unknown sour ce ellipticity and 
the o bservational uncertainties (See e.g. ISchneider et al 
2000). In addition, we also measured the ^-statistic and 
£ c -statistic densitometry: 

C(r; r max ) = R{r' < r) - R(r < r' < r max ) 

2 /"''max 

(7t(r))dlnr, (19) 



1 - (rmaxA) 2 

Cc(r; r 2 ,r max ) = R(r' < r) - R(r 2 < r' < r max ) 
= [l - (r/r max ) 2 ] ((r;r max ) + 
1 - {r 2 /r max ) 2 



(^max/^) 2 



1 



C(r 2 ;r max ) 



(20) 



where (7t(r)) is the averaged tangential ellipticity inside 
an annulus (7rr 2 E cr i t K(r) is the mass enclosed within the 
radius r). In contrast with the £ statistics, ( c can be di- 
rectly compared with R since they only di ffer by a con- 
stant value that does not change with radius l) Clowe et all 
|2000|) . We used r 2 = SUh^kpc and r max = 744/i" 1 kpc. 
The scaling factor for the mass has been derived from the 
UBVRIJK photometric redshifts of sources. Background 
galaxies have been selected in the magnitude range I < 24 
and cluster galaxies have been rejected using a photo- z se- 
lection. Moreover, we considered background galaxies with 
z p h > 0.4 for which the lensing signal is significant. The 
limiting magnitude was chosen in order to compromise 
between the depth, which defines the galaxy number den- 
sity, and the need for a good estimate of the source red- 
shift distribution. Since our source population is similar to 
IVan Waerbeke et alJ l|2002t) . we checked our redshift his- 
togram has the same shape 10 as their sample. Both sam- 
ples turned out to be similar, so we finally used their pa- 
rameterized redshift distribution, because it is based on a 
larger sample than ours. With this requirement, the weak 
lensing signal directly makes a test on the reliability of 
strong lensing models extrapolations beyond the Einstein 
radius. 



after subtraction of the cluster population 
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Table 4. Single halo models. The Einstein radius is the same in the best models for strong lensing : Re — 50 kpc. 
For our two strong lensing models, the total mass inside this radius is M(r < Re) — (2.8 ±0.1) x 10 13 /i _1 M Q . Errors 
du e to uncertain t ies in arcs redsh ift are omitted . Also reported p revious works results fo r compar ison. Me93 refers 
to llMellier et allll993l). Ha97 to fHammer et all Il997[) . Mi95 to l|Miralda-Escudel Il995[) EF99 to jEttori fc Fabiari 
and A102 to l|Allen et all l2002t) . When known, the authors' values are recomputed in our adopted cosmology 
and with the 1.6 sources redshifts. For both papers, the center of potential location is assumed to match the center of 
cD or is not reported. Me93 core radii have been scaled in order to take into account the departs between their profile 
and an exact softened isothermal sphere. As well, Ha97 find a slope (3 w 0.85 instead we have only considered models 
with (3 = 1 (see Eq. (J2J)). A102 uses a NFW profile and only gives the scale radius value but we report on the same line 
our own measured values for ellipticity and position angle from Chandra X-ray brightness. The third column r c \r s 
corresponds either to the scale radius either to the core radius. Here, we convert all the position angles in a common 
definition, which is clockwise from North to East. The original paper do not report angles in the same frame but we 
made the correction except for Ha97 for which we do not know what is the reference. But in any case, the position 
angle is so constrained that these authors must have found a similar orientation as the other ones. Our definition is 
more valuable and self-consistent between Chandra, ROSAT, VLT and HST data. Models labels with a S refer to 
purely strong lensing modeling whereas a W stands for purely weak lensing fits. The last row (cD+DM) concerns 
the last family of profiles with a cD and dark matter halo components and which is simultaneously constrained by 
strong+weak lensing. In the first column, we report the permitted inner slope for generalized NFW profiles (see Sect. 



model 


(km/s) 


( h" 1 kpc) 




>"200 

( h' 1 kpc) 
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0.25 ± 0.05 


59 ± 1 


0.2 ± 0.4 


0.2 ± 0.4 


W-isoT 
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< 45 


















Me93 


~ 1000 


4.5 - 7 










0.15 - 0.33 
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Mi95 


~ 1200 


~ 10 










~ 0.22 


~ 58 






Ha97 


~ 1100 


5-10 










0.18 - 0.23 


?47.5 ± 5? 






EF99 


~ 930 


~ 90 




r&oo ~ 540 














A102 




107 - 120 










~ 0.20 


58 ± 7 






cD+DM 


0.8 < a < 1.1 


~ 85 


~ 0.6 








0.22 ± 0.06 


58 ± 2 


0! 


0! 



Fig. El shows the radial mass profile of the best IS 
and NFW models. The projected mass density has been 
averaged inside circular annuli. As expected, the two 
best fits are quite similar between the two critical radii. 
Discrepancies only appear in the innermost and outermost 
regions. However, the shear profile derived from the VLT 
data fails to disentangle the models built from strong lens- 
ing. Both are consistent with the signal down to the virial 
radius r 2 oo ~ l/i. _1 Mpc. Tabled lists the values of the 
best fit parameter set for the weak lensing analysis. It 
is in good agreement with the inner strong lensing mod- 
els, though the total encircled mass is smaller. The con- 
straints on the concentration parameter are weak and a 
broad range of values are permitted. However, a low value 
similar as expectations for clusters is still marginal and 
surprisingly the weak lensing analysis also converges to- 



ward a rather larger concentration. This discrepancy with 
cluster expectation values, even when using together weak 
and strong lensing constraints, shows that the global prop- 
erties of the potential well are hard to reconcile with a 
simple NFW mass profile. However, if the contribution of 
the cD stellar mass profile strongly modifies the innermost 
mass distribution of the cluster and significantly contami- 
nates the concentration parameter inward, our statement 
based on strong and weak lensing models might be wrong. 
We therefore single out the cD potential and add its con- 
tribution to the model and we included the fifth image 
parameters in order to probe the very center where the 
cD mass profile might play an important role. 
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Fig. 5. Projected mean surface density. Solid: NFW best fit, dashed: IS. Notice how models match near the arcs 
locations. Between 2.5 and 28 arc-sec, the mean convergences R —< k > (0, r) differ from each other by less than 3 
percents. Faint discrepancies appear in the outer regions through weak lensing signal and at the very center. Constraints 
due to the fifth image candidate are detailed in Sect. 13. 431 The weak lensing data are deduced from the £ c -statistic. 



3.4. The cD+DM mass profiles constrained with the 
fifth image. 

We now consider a two-component mass profile : an 
inner stellar component attached to the cD galaxy and a 
cluster dark matter halo. The fifth central image will con- 
tribute to constrain the innermost lens model, whereas the 
external arcs and the weak lensing profile should constrain 
most of the outer cluster halo. 

3.4.1. Centering the lens with the fifth image 

Before introducing the stellar component, let us check 
the influence of the fifth image knowledge on the centering 
of a single DM potential. Fig. shows the permitted area 
for the DM potential center relative to the cD. The con- 
tours on the top are the expectations for the IS and NFW 
models, if the fifth image is not taken into account. The 
offset with respect to the cD centroid is 0.22" West, but 
the contour ellipses are of size 1.1" x 1.6". Nevertheless, 
the assumption that the center of cD galaxy coincides with 
the cluster center is consistent with the data. When the 
fifth image is added (mainly its brightest knot), the con- 
tours shrink by a factor of 2 in size, as shown in Fig Eh, 
but still keep the central cD position inside, with a small 
offset with respect to the cD light centroid of 0.16" West 
for the IS model, and 0.18" West for the NFW profile. The 
box sizes of permitted positions are much smaller ellipses 
of about 0.6" x 0.9". Since these error boxes are about the 
size of the uncertainties of the cD centroid position (see 
Tabic in the following we will then assume the cD is 
centered on the cluster center. It is worth noticing that 
even with the significant reduction of error bars provided 
by the fifth image, the residual uncertainty on the cen- 
troid position of the lens may in principle permit to both 



IS and NFW to fit the lensing data if we do not assume 
the cD center is not exactly on the cluster center. 

3.4.2. Modeling together the stellar and DM mass 
profiles 

The properties of the lens configuration (including the 
fifth image) provide enough constraints to attempt a mod- 
eling which will probe clear differences between observa- 
tions and IS/NFW predictions. The deflection and the 
magnification of the NFW model are smaller than for an 
IS one. We expect the fifth image to show a difference of 
0.2" in position and 0.75 in magnitude. The observations 
and IS/NFW predictions reported in table El and Fig. 0] 
already show a trend which supports a flat-core model 
against a cuspy NFW profile. 

The following analysis uses together the fifth image prop- 
erties, the weak lensing data and the giant arcs in order 
to constrain the shape of the innermost mass profile. We 
also add the cD stellar contribution to the overall mass 
because it is no longer negligible at the very center. Note 
that the generalized NFW models expressed in Eq. has 
a free parameter a. Its projection is reported in Eq. (|B.6(I . 
In more details: 

— The fifth image central coordinates are introduced be- 
cause of their constraints on the central lens modeling. 
In fact, the brightest knot in the large arc AO is re- 
quired to correspond to the brightest detected spot at 
the center r$ ss 0.9. 

— The center of the cD is precisely the center of the po- 
tential well. 

— The dark matter halo is modeled by an elliptical mass 
distribution: a softened IS profile or a generalized- 
NFW profile, as expressed in Eq. (|B.6|I . Hence, the 
IS profile has 4 free parameters, namely r c , a v , e, PA 
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Fig. 6. Permitted potential center region (Ax 2 /N = 
1,2,3) relative to the cD center (0,0) coordinates. Solid 
(resp. dashed) contours refer to the NFW (resp. IS) model. 
The upper panel displays contours without the fifth image 
location knowledge and the lower using the 5th image can- 
didate discussed in Sect. 12.51 as an additional constraint. 
The crosses (resp. diamonds) correspond to the \ 2 mini- 
mum location for the NFW (resp. isothermal) modeling. 
Note the great enhancement induced and the significant 
shift in the NFW case compared to the upper panel. The 
deflection angle is merely too small to be consistent with 
the fifth image. 

whereas the generalized- NFW models has 5 : a, r s , k s , 
e, and PA. 

— We model the stellar component with an Hernquist 
profile (p(x) = x _1 (l + a;) -3 ) of the projected form : 



1 - F(y) 

v 2 -l 



(21) 



with F defined in Eq. I|17|). y = r/r g and r g sa 
7.2 h^ 1 kpc a scale radius, related to the I band 

luminosity through the relation: 



\ =0.1 IT. 



27TS 2 S cr jf \ L 



where T = M/Lj. The stellar component is elliptical 
and has the central ellipticity e* = 0.15 and orientation 
PA* = 69° deduced from the light distribution in the 



I band. The cD stellar mass-to-light ratio T is the last 
new free parameter of the model. 

— The weak lensing x 2 term (Eq. IT%|l is directly added 
to the purely strong lensing \ 2 defined in Eq. (JJJ). 

3.4.3. The inner slope of the DM halo 

The constraints provided by the fifth image can 
roughly be illustrated as follows. If r 5 w 0.9 corresponds 
to a source position u w 1.85" 11 deduced from the sin- 
gle component outer arcs. The lens equation reads: u = 
7-15 1 1 — K,(r 5 )\ with the bending angle R{r) = ~-^jr. Hence, 
the averaged convergence within the fifth image candidate 
radius which is plotted on Fig [5] is 

2.9 < k(< r 5 ) < 3.7. 

It is now possible to test the overall permitted central 
contribution by mapping the global \ 2 m the r c — M/Li 
or a — M/Lj spaces, after marginalization over the other 
parameters. For a reasonable mass-to-light ratio w 2 — 3, 
one can see on Fig that : 

— Both IS and generalized-NFW models rule out M/Lj 
of the cD stellar component larger than 9 at a 3-cr 
confidence level. Its value preferentially ranges within 
1 to 5. 

— The softened IS profile still provides the best model. 
It is also consistent with the strong lensing data with 
few variations of parameters compared to the single 
component modeling. The introduction of the stellar 
component does not introduce large variations in the 
best fit parameters set compared to the single dark 
halo modeling. 

— The cuspy models have a narrow permitted range of 
slope which is centered around a = 1 of the NFW 
model. Note that while the position of the fifth image 
provides interesting boundaries on the cluster center 
position, it does not provide constraints on the slope 
a. Including the fifth image only reduced the a upper 
limit by 10% and does not changes its lower bound. For 
reasonable values of M/Lj, we find 0.8 < a < 1.1 (2- 
<t) . This range excludes very low values of a and seems 
to contradict the fact that IS with flat core better fit 
the data than generalized-NFW models over the whole 
l/?, _1 kpc < r < 600 ft. -1 kpc range, even for a large 
amount of stellar mass. 

4. Discussion 

4.1. The radial mass profile of MS2137. 3-2353 

The exceptional data set allowed us to constrain the 
density profile over three orders of radius ranging from 2 
to 700 hr 1 kpc. Despite the fact that weak lensing data do 
not cover a wide enough range in order to reveal its full 

11 Few variations are observed when modeling with the NFW 
or the IS model. We also neglect ellipticity terms which have 
a weak importance near the center 
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straining second order moments of the fifth image probing 
both convergence and shear inside 1 arc-sec radius. 

Nevertheless, the new constraints provide three levels 
of information concerning the MS2137. 3-2353 radial mass 
profile. 

(i) The simple use of the radial and tangential arcs 
systems without the fifth image neither the outer 
shear data cannot disentangle between NFW-like and 
isothermal profiles. A large family of cuspy and flat 
models are consistent with these data. (See Figs.|Hland 
Eland appendix lB|) 

(ii) The combined weak lensing and arcs (i) data tell us 
that either isothermal profiles either NFW-like profile 
with 0.8 < a < 1.1 are permitted. 

(iii) Actually, the fifth image knowledge favors flat cores 
(xjs ~ 3-8) but puts strong limits of the couple a — T 
for NFW-like models (xLi -5.1). 



Fig. 7. Two components x 2 contours. Upper panel : con- 
straints on the couple r c — M/Lj for the IS model. Note 
the small modifications due to the introduction of the stel- 
lar component. Lower panel : constraints on the couple 
a — M/Lj for the pseudo-NFW family profiles. One can 
see an overlay of the same confidence regions when the 
fifth image position is known (solid) and when it is not 
(dashed). The inner slope a ranges between 0.8 and 1.1 
for a reasonable mass of stars (T/ < 5). The lower panel 
rejects very shallow profiles (a < 0.6) whereas one can see 
on the upper one that flat cores provide good fits. 



efficiency 12 , we performed a self consistent modeling of 
the critical strong and sub-critical weak parts of the lens- 
ing cluster MS2137. 3-2353. At the other side, it is worth 
noticing that the improvement provided by the fifth image 
is still under-exploited because of the poor resolution of 
its shape. The location of its brightest spot only provides 
constraints on the center position of the lens and on the 
overall enclosed mass (by the way, revealing a degeneracy 
between T and a). A good knowledge of the magnification 
and shear would be able to break this degeneracy by con- 



All together, the new constraints are in good agree- 
ment with isothermal model with flat core and rule out 
generalize d-NFW models with slopes as steep as thos e pro- 
posed by l|Moore et all Il998t iGhigna et all l2000j) . The 
slope range found for the generalized-NFW profiles can 
be easily explained and correspond to expectations. The 
calculations detailed in appendix iBl fFias. IS1 and show 
how the knowledge of the lensing configuration, as de- 
rived from giant arcs and the shear field, bounds the free 
parameters for cuspy-NFW profiles. The knowledge of the 
critical lines radii, the weak lensing at intermediate scales 
as well as the length of the radial arc are introduced in 
order to fix semi-analytically r s , a and k s . 




12 the ("-statistic S/N ratio at radius r goes like In r max /r for 
an isothermal profile. 



Fig. 8. Plot of various ( c curves taken along the relation plot- 
ted in blue (light gray) on FigQUthat give the observed critical 
lines. Also plotted the observed Cc( r ) deduced from weak lens- 
ing. Only a small scatter around a = 1 is allowed. We also 
report the best fit isothermal profile (dashed). 

A lower curvature for NFW-like profiles can explain 
the apparent paradox discussed at the end of Sect. 13. 4131 
Flat softened isothermal profiles are favored whereas low 
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values of a are rulcd-out. This trend was reported by 
iMiralda-Escudel l|l995|) when he tried to fit the dark mat- 
ter halo with a density profile of the form 

p(x) oc af x (l + x)~ 3/2 or p(x) (x x^ 1 (1 + xy 1 . 

Hence, a simple smooth modification from the inner slope 
a <C 1 to an asymptotic p oc r -3 outward cannot match 
the whole lensing data. Sharper changes must occur at 
a small radius which behaves as an effective core radius, 
leading to a high curvature close to the radial arc. The 
scale radii derived from the marginalization of Fig [7| are 
very small and still scatter around the r s w 90ft. _1 kpc 
value obtained from the single component NFW. The sen- 
sitivity of models to radial curvature is clearly visible on 
Fig El 




a? a-1 c.e i" fi ■ o 



Fig. 9. Relations in the parameter space a — r s deduced from 
the critical lines Eq. lB.4at and Eq. <B,4bll (blue or light gray 
region), from the ("—statistic and Eq. (|B.4a(l (red or dark gray) 
and from the radial arc length relative to its counter-image 
(dashed yellow area). Errors take into account variations of el- 
lipticity, of the mass-to-light ratio and observational uncertain- 
ties. The resulting permitted region is well consistent with what 
is found from the modeling. One can see that 0.7 < a < 1.1. In 
other words, the halo density profile must be much shallower 
than the value 1.5 suggested by recent numerical simulations 
jMoore et allll99&lGhigna et al"Ll2000l) . The permitted range 
for r s is restricted to small scale radii r s ~ 100 h~ kpc. 

The whole best fit generalized-NFW profiles show a 
high concentration for the dark matter halo. This trend 
is confirmed by weak lensing up to 700 /i -1 kpc, in con- 
trast with other weak lensing cluster analyses which find 
smaller concentration than ours, but more consistent with 
numerical CDM simulations. The role of stars does not 
change this conclusion. So, if generalized-NFW models 
are acceptable, it is important to confirm that in the 
case of MS2137. 3-2353 they imply the concentration to 
be stronger than numerical predictions, ft is therefore 
important to confirm these results by using a different 



method. Recently, ISand et all l)2002|) have reported com- 
parable slope constraints using the velocity dispersion of 
stars at the center of the cD and the positions of crit- 
ical lines. Conversely, any lensing model should be con- 
sistent with the information on the kinetic of stars they 
measured, so it is necessary to compare our predictions 
with their data. Nonetheless we plan to show elsewhere 
(Gavazzi et al., in preparation) that the velocity disper- 
sion usually measured from the FWHM of absorption lines 
in the galaxy spectrum no longer hold if the distribution 
function of stars is far from a Maxwellian as mentioned in 
iMiralda-Escudl l|l995|h 

Finally, we checked that the introduction of galaxy 
halo perturbations under the form of massive haloes at- 
tached to the surrounding cluster galaxies does not change 
our conclusions. Such perturbations have poor conse- 
quences for the weak lensing results but are likely affect 
slightly the fifth image location. We show in appendix [O 
that a significant modification of the fifth image due to 
galaxy halo perturbations implies to put a huge mass on 
each galaxy. Such an amount of mass would destroy the 
quality of the arcs fit. 

4.2. Effects of non constant ellipticity and isodensity 
twist on the radial arc 

At the tangential arc radius (r ~ 50 /i -1 kpc) we 
measure a robust offset angle A6 = 13° ± 4° between 
the diffuse stellar component and the DM potential 
orientation. This result is confirmed by the Chandra 
X-rays isophotes contours as shown on Fig [3 Previous 
strong lensing modelings in the presence of important cD 
galaxies never clearly established such a behavior because 
the uncertainties of the models obtained with tangential 
arcs only were too loose for the isopotential orientation. 
However, f or the nearby elliptical galaxy NGC720 
iBuote etaH l)2002i) and iRomanowskv & KochanelJ 
(1993) 13 studied such a misalignment between the light 
distribution and the surrounding dark halo revealed by 
X-ray emissivity. RK showed that the stellar misalign- 
ment can be explained by a projection effect of triaxial 
distributions with aligned main axis but different axis 
ratios. 

Moreover all the best fit modelings show a tiny but 
robust remaining azimuthal offset (< 0.3") between the 
modeled radial arcs (Al and AT) and the position actu- 
ally observed on the HST image (highlighted on recon- 
structions of Fig^l. 

We verified that it is not due to a bad estimation of the 
source position since any small source displacement pro- 
duces a large mismatch between the counter-arc A5 and 
the model. This pure azimuthal offset led us to investigate 
the possible effect of a variation of the ellipticity and po- 
sition angle of the projected potential close to the radial 
arc radius. Such a trend is also favored by an increase of 
ellipticity on the X-ray isophotes with radius. 

13 see references therein for a list of analogous objects 
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If we had implemented the availability of using models 
with a variation of orientation and cllipticity with radius 
in the inversion software, we would have found that the 
orientation of the potential major axis tends to the ori- 
entation of stars (see Fig El) when looking further in. At 
the same time the potential becomes rounder. We roughly 
checked this behavior by modeling the lens configuration 
with two distinct (and discontinuous) concentric areas (in- 
side and above 8 arc-sec). The rays coming from the source 
plane and giving rise to the outer arcs AO, A2, A4, A5 do 
not suffer any inner variation of the potential symmetry 
(provided that the overall mass inclosed in the Einstein 
radius remains the same). Hence the previous modelings 
remain valid for the outer parts whereas the inner can be 
twisted and made rounder in order to alleviate the off- 
set problem. A small twist A PA ~ 8° in the direction of 
stars gone with a smaller ellipticity (e ~ 0.2) turns out 
to suppress efficiently the azimuthal offset near the center 
without affecting the external arcs. This analysis is not 
exhaustive in the sense that maybe different explanations 
can be found but its main virtue is to show that high spa- 
tial resolution like HST imaging of numerous multiple arcs 
makes a lens modeling so binding that it becomes possible 
to extract much more information than the simple fit of 
elliptical models. In addition with the hot ICM properties 
(see e.g. iRomanowskv fc Kochanekl Il998h . we could cer- 
tainly start more detailed studies of potential with twist 
effects and eventually start to probe the triaxiality of dark 
matter halos if we can observe a large number of multiple 
arc systems in clusters. 

These results strengthens the argument of lMiralda-Escudel 
(2002) upon which the ellipticity of DM halos makes in- 
consistent the hypothesis of self-interacting dark matter. 

5. Conclusion 

By using strong and weak lensing analysis of HST and 
new VLT data of MS2137. 3-2353 we found important new 
features on the lensing configuration: 

The photometric redshifts or the radial and the tan- 
gential arcs are both at 1.6 ± 0.1 in excelle nt agreement 
with the recent spectroscopic observations of ISand et all 
The extraction of the cD diffuse stellar light has permit- 
ted to detect only one single object which turns out to be 
at the expected position of the fifth image. Furthermore, 
its orientation, its ellipticity, its signal-to-noise ratio and 
its morphology correspond to those expected by the lens 
modeling. Unfortunately, the poor determination of its 
shape properties hampers the use of its geometry as a local 
estimate of the magnification matrix toward the center. 
Using the fifth image together with the weak lensing anal- 
ysis of VLT data, we then improved significantly the lens 
modeling. The radial mass profile can then be probed over 
three orders in distance. This additional constraint seems 
to favor isothermal profiles with flat core or generalized- 
NFW profiles with 0.8 < a < 1.2 without introducing 
the fifth image knowledge. When this constraint is added 
together with the prior motivated that cD center and clus- 



ter halo center are the same, we favor flat core softened 
isothermal spheres. The position of the fifth image is in 
better agreement with an isothermal model than an NFW 
mass profile. In addition, it is worth noticing that the ki- 
netics of stars should be analyzed in details, considering 
a precise distribution function that may depart from the 
commonly assumed Maxwellian. 

We point out a misalignment between the diffuse stel- 
lar component major axis and both the lens potential and 
the X-ray isophotes. We argue it is produced by the triax- 
ial shape of the mass components. This extends the previ- 
ous demonstration of the cllipticity of the projected dark 
matter halo. This work is a first attempt to improve strong 
lensing observables and modeling in order to probe both 
the central DM cusp/core and the triaxiality of DM halos. 

It would be essential to confirm the detection of the 
fifth image. Fig. [3| shows that the spectral energy distri- 
bution expected for the fifth image is different than the old 
stars dominated cD emission. We therefore expect the fifth 
image to show up on an optimal image subtraction U — XJ 
(A being optimized). We attempted to use this technique 
on our present data but the poor resolution (~ 0.6") on 
the U and J ground based images prevent any significant 
enough detection. We conclude that only a high resolution 
observation with the Space Telescope in UV-blue wave- 
lengths or in a peculiar emission-line is among the best 
constraints one could envision in the future. 

There is not yet evidence that similar studies as this 
work can be carried out on other ideal lens configurations. 
The strength of the diagnostic on the radial mass profile is 
however so critical that we must apply this technique to a 
large sample in order to challenge collision-less CDM pre- 
dictions on a realistic number of clusters of galaxies with 
eventually a test of the role of dominant central cD galax- 
ies. The simultaneous use of weak lensing data should be 
more relevant for wider fields in order to check also a r -3 
fall-off on the density profile predicted at large distance 
by CDM simulations. 

Acknowledgements. The authors would like to thank Jordi 
Miralda-Escude for numerous and fruitful discussions and com- 
ments on the content and the outlooks of that work. We also 
thank Peter Schneider, Tom Broadhurst and Avishay Dekel 
for useful comments. Special thanks to Stella Seitz who kindly 
provided their B and R VLT images of the cluster and for 
her comments on that work and to Chuck Keeton who made 
the last version of the lensmodel software available. The au- 
thors also thank the referee for his comments that help to clar- 
ify several points of the paper. This work was supported by 
the TMR Network "Gravitational Lensing: New Constraints 
on Cosmology and the Distribution of Dark Matter" of the EC 
under contract No. ERBFMRX-CT97-0172. 

References 

Allen, S. W. 1998, MNRAS, 296, 392 

Allen, S. W., Schmidt, R. W., & Fabian, A. C. 2002, MNRAS, 
334, Lll 

Arabadjis, J. S., Bautz, M. W., & Garmire, G. P. 2002, ApJ, 
572, 66 



16 



R. Gavazzi et al.: Radial mass profile of MS2137. 3-2353 



Athreya, R. M., Mellier, Y., van Waerbeke, L., et al. 2002, 

A&A, 384, 743 
Bartelmann, M. 1996, A&A, 313, 697 
Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 
Blumenthal, G. R., Faber, S. M., Flores, R., & Primack, J. R. 

1986, ApJ, 301, 27 
Bolzonella, M., Miralles, J.-M., & Pello, R. 2000, A&A, 363, 

476 

Bode, P., Ostriker, J.-P., Turok, N. 2000, ApJ, 556, 93 
Bruzual, A., G. & Chariot, S. 1993, ApJ, 405, 538 
Bullock, J. S., Kravtsov, A. V., & Weinberg, D. H. 2001, ApJ, 
548, 33 

Bullock, J. S., Kolatt, T. S., Sigad, Y., et al. 2001, MNRAS, 
321, 559 

Buote, D. A., Jeltema, T. E., Canizares, C. R., & Garmire, 

G. P. 2002, ApJ, 577, 183 
Burke, W. L. 1981, ApJ, 244, LI 

Clowe, D., Luppino, G.A., Kaiser, N., Gioia, M.I. 2000, ApJ, 
539, 540 

Clowe, D., Schneider, P. 2001, A&A, 379, 384 
Dalai, N. & Kochanek, C. S. 2002, ApJ, 572, 25 
Ettori, S. & Fabian, A. C. 1999, MNRAS, 305, 834 
Fort, B., Le Fevre, O., Hammer, F., & Cailloux, M. 1992, ApJ, 
399, L125 

Gavazzi, R. 2002, New Astronomy Review, 46, 783 
Ghigna, S., Moore, B., Governato, F., et al. 2000, ApJ, 544, 
616 

Gioia, I. M., Maccacaro, T., Schild, R. E., et al. 1990, ApJS, 
72, 567. 

Haehnelt, M. G. fc Kauffm ann, G. 2002, submitted to 

MNRAS., astro-ph/0208215 
Hoekstra, H., Franx, M., Kuijken, K., & van Dokkum, P. G. 

2002, MNRAS, 333, 911 
Hammer, F., Gioia, I. M., Shaya, E. J., et al. 1997, ApJ, 491, 

477 

Jing, Y. P. & Suto, Y. 2002, ApJ, 574, 538 

Keeton, C. R. 2001, ApJ, 561^6 

Keeton, C. R. 2001, preprint, astro-ph/0102340 
Keeton, C. R. 2001, preprint, astro-ph/0111595 
Kelson, D. D., Zabludoff, A. I., Williams, K. A., et al. 2002, 
ApJ, 576, 720 

Klypin, A, Lravtsov, A, Valenzuela, O. & Prada, F. 1999, ApJ, 
523, 32 

Kneib, J. P., Mellier, Y., Fort, B., & Mathez, G. 1993, A&A, 
273, 367 

Kneib, J.-P., Ellis, R. S., Smail, I., Couch, W. J., & Sharpies, 

R. M. 1996, ApJ, 471, 643 
Mailer, A. H. & Dekel, A. 2002, MNRAS, 335, 487 
Mellier, Y. ARAA 37, 127. 

Mellier, Y., Fort, B., & Kneib, J.-P. 1993, ApJ, 407, 33 

Merritt, D. 1985, AJ, 90, 102 

Merritt, D. 1985, MNRAS, 214, 25 

Metcalf, R. B. & Madau, P. 2001, ApJ, 563, 9 

Milosavljevic, M. , Merritt, D., Rest, A., & van den Bosch, 

F. C. 2002, MNRAS, 331, L51 
Miralda-Escude, J. 1995, ApJ, 438, 514 
Miralda-Escude, J. 2002, ApJ, 564, 60 

Moore, B., Governato, F., Quinn, T., Stadel, J. & Lake, G. 

1998, ApJ, 499, L5 
Moore, B., Ghigna, S., Governato, F., et al. 1999, ApJ, 524, 

L19 

Moore, B., Gelato, A. J., P earce, F. R., Quilis, V. 2001, 
preprint, astro-ph/0002308 



Natarajan, P., Kneib, J., & Smail, I. 2002, ApJ, 580, Lll 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 
493 

Navarro, J. F. & Steinmetz, M. 2000, ApJ, 528, 607 

Osipkov, L. P. 1979, Soviet Astronomy Letters, 5, 42 

Pello, R., Bolzonella, M., Campusano, L. E., et al. 2001, 

Astrophysics and Space Science Supplement, 277, 547 
Romanowsky, A. J. & Kochanek, C. S. 1998, ApJ, 493, 641 
Salucci, P., Burkert, A. 2000, ApJ, 537, L9 
Sand, D. J., Treu, T., & Ellis, R. S. 2002, ApJ, 574, L129 
Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational 

Lenses, XIV. Springer- Verlag. Also Astronomy and 

Astrophysics Library. 
Schneider, P., King, L., & Erben, T. 2000, A&A, 353, 41 
Smith, G. P., Kneib, J.-P, Ebeling, H., Czoske, O., & Smail, 

I. 2001, ApJ, 552, 493 
Sofue, Y. & Rubin, V. 2001, ARA&A, 39, 137 
Stocke, J. T., Morris, S. L., Gioia, I. M., et al. 1991, ApJS, 76, 

813 

Spergel, D. N., & Steinhardt, P. J., 2000, PRL, 84, 3760. 
Tyson, J. A., Kochanski, G. P., Dell'Antonio, I. P. 1998, ApJ, 
498, L107 

Van Waerbeke, L., Mellier, Y., Pello, R., et al. 2002, A&A, 393, 
369 

V erde, L., J imenez, R. 2002, submitted to MNRAS, 

astro-ph/0202283 
Willick, J. A. fc Padmanabhan, N. 2000, preprint, 

|astro-ph/0012253| 
Wise, M. W. & McNamara, B. R. 2001, Two Years of 

Science with Chandra, Abstract from the Symposium held 

in Washington DC, September 2001, 217 



Appendix A: Table of identified dots. 

Appendix B: Analytical constraints on 
generalized- NFW profiles 

We consider the simplified case where the lens is de- 
scribed by an elliptical density profile which has a small 
ellipticity e and a small stellar contribution. We neglect 
terms in the multi-polar development higher than the 
quadrupole (first order ine). 

One can easily write a set of equations that the system 
must verify : the critical lines locations, the lens equation 
relating the radial arc Al and its counter-image A5. One 
can also force the model to fit the weak lensing constraints 
at large radius, say r w = 160 ft. -1 kpc. The tangential line 
is known to pass by the point 

r t = [-11.25"; 10.25"] -> r t = 15.2", ip t = 137.7°, 

the radial line to pass by the point 

r r = [-3.16"; 3.93"] -> r r = 5.0", tp r = 137.6°, 

the associated point in A5 is 

r c = [8.3"; -22.4"] -> r c = 23.9", tp c = -69.5°, 

and the weak lensing (— statistic constraint reads 

n(r w ) = 0.1 ± 0.03. (B.l) 
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Table A.l. X and Y coordinates in arc-sec for the 26 knots used in the modeling. Coordinates are referred to the 
cD center and are oriented as in FigQ] The first part consists of the tangential arc system. The first 13 are 5 times 
imaged and the next 6 are triply imaged. The (*) denotes the brightest spot in the arcs and is the only point which 
are seriously observed in the fifth image. On the other side, the radial arc consists on triply images dots only. Alin 
(resp. Alout) is the inner (resp. outer) part of the composite radial arc. Both are imaged into the Eastern extremity of 
A5. In this later case, associations of mid points (x,y)=(8.3,-22.4) are somewhat ambiguous and are given less weight 
for the modeling. The tangential and radial systems were used simultaneously since we established they are at the 
same redshift z s ,phot — 1-6. The system {AT-A6} is found at a slightly lower photometric redshift and can only be 
compared with the modeling at a later stage. 



ID 


A02 


A01 


A2 


A4 


5th 


1 


-14.18 


6.39 


-10.64 


11.67 


12.06 


5.24 


-2.65 


-18.99 


0.63 


0.70 


2 


-14.26 


5.48 


-8.65 


13.24 


12.29 


5.22 


-2.49 


-18.79 


0.65 


0.77 


3 


-14.28 


5.74 


-9.62 


12.58 


12.14 


5.38 


-2.38 


-18.94 


0.60 


0.67 


4 


-13.51 


7.96 


-11.39 


10.92 


11.98 


5.48 


-2.46 


-19.13 


0.61 


0.70 


5 


-13.95 


5.89 


-8.35 


13.27 


12.11 


5.76 


-1.81 


-18.99 


0.61 


0.63 


6 


-14.00 


5.22 


-5.81 


14.54 


12.33 


6.12 


-0.78 


-19.25 


0.60 


0.60 


7 


-13.74 


5.55 


-5.32 


14.63 


12.21 


6.56 


-0.26 


-19.20 


0.54 


0.61 


8 


-13.71 


5.39 


-5.10 


14.80 


12.38 


6.45 


-0.08 


-19.14 


0.50 


0.61 


9 


-13.58 


6.05 


-5.73 


14.29 


12.01 


6.82 


-0.23 


-19.32 


0.51 


0.65 


10 


-12.65 


7.78 


-6.98 


13.45 


11.52 


7.40 


-0.24 


-19.57 


(*)0.59 


0.66 


11 


-12.43 


8.16 


-7.24 


13.26 


11.50 


7.22 


-0.71 


-19.56 


0.57 


0.66 


12 


-11.89 


9.21 


-9.08 


12.07 


11.30 


7.22 


-1.07 


-19.57 


0.60 


0.73 


13 


-11.13 


10.42 


-11.60 


9.91 


11.48 


6.99 


-1.24 


-19.42 


0.62 


0.76 


14 










11.14 


7.06 


-1.48 


-19.61 


0.62 


0.65 


15 










11.12 


6.84 


-1.88 


-19.63 


0.65 


0.67 


16 










11.27 


6.02 


-3.10 


-19.36 


0.70 


0.66 


17 










10.69 


5.90 


-3.92 


-19.63 


0.79 


0.79 


18 










10.83 


5.62 


-4.31 


-19.06 


0.73 


0.84 


19 










11.14 


5.07 


-4.76 


-19.41 


0.77 


0.82 



ID 


A5 


Alin 


Alout 


20 


7.9 


-22.4 


-1.8 


3.0 


-4.1 


5.4 


21 


8.3 


-22.4 


-2.1 


3.2 


-3.9 


5.1 


22 


8.3 


-22.4 


-2.3 


3.4 


-3.6 


4.6 


23 


8.3 


-22.4 


-2.5 


3.5 


-3.6 


4.3 


24 


8.3 


-22.4 


-2.7 


3.6 


-3.2 


4.1 


25 


8.5 


-22.3 


-2.8 


3.8 


-3.0 


4.0 


26 


10.1 


-21.4 


-3.2 


3.2 


-3.8 


3.8 



In the following, subscripts t, r and c denote the values 
taken at position r t , r r , r c respectively. If we write the 
magnification matrix as: 



1 - K + 7i -72 
-72 1 - K - 7i 



(B.2) 



and the lens equation between radial arc Al and its 
counter-image A5 position: 



u = r r - V(f){r r ) =r c - V<^(r c ) 



(B.3) 



Reduced to the first terms in e, the equation of the 
tangential (k + 71 = 1), the radial (n — 71 = 1) lines and 
the radial relation between Al and A5 respectively yield: 

«t + [36 - Kt}e t = 1-St (B.4a) 
2n r — R r — [3£ r — K r + r r K (r r )]e r = 1 — S r (B.4b) 

k c + r\R r + e c (£ c - «c) + e r ?7(£r — K r ) = I + r) — S c , 

(B.4c) 



where e l=t , r ,c = ecos(2(ipi - tp )), k' = rj = r r /r c 
0.21 and 



£(r) = 



dr' r' 3 K(r') 



(B.5) 



(fully detailed in iMiralda-Escudet Il995ft . The terms 
Si=t,r,c are the small corrections from the stellar contri- 
bution. St ~ K«(r t ), S r « 2re*(rv) — «*(r r ) and S c w 
rjR^{r r ) + K*(r c ) which are of a few percents order and 
scale like the mass-to-light ratio T. 



If we now project a general 3D density profile (QJ into: 



k(t) =2n s x 1 - a {{l + x) a - 3 + 
»i 



(3-a) / dy(y + x) a -*(l- Vl^f)} 



R{r) = An s x^ a {^- 2 Fi(3 - a, 3 - a, 4 - a; -x) 



(B.6) 



+ C Ay (y + x) a - 3 - ^ l y2 } 
Jo V 

with x = r/r s and k s = p s r s /E cr i t we can constrain all 
the parameters r s , k s and a for a given ellipticity and a 
given mass-to-light ratio T. We retrieve the NFW profiles 
for a = 1. We also need to assume a position angle and 
an ellipticity that we set equal to the values deduced from 
the modeling: ipo — 5° , e = 0.24. We analyzed departs 
from this value. 

In fact, we solve the set of Eq. l|B.4a|) and Eq. l|B.4b|l 
for r s and function of the inner slope a. Notice 

that the whole set of Eq. (|B.4fl would in principle be 
sufficient for constraining exactly the triplet [a, r s ,K s ]. 
Nevertheless, the radii inferred in these equations are very 
similar and thus the solution suffers a high sensitivity to 
the uncertainties on the values of r r , r c and r t . 

The numerical modeling deals with much more con- 
straints than the relations Eq. HB.4(I and Eq. IjB.ip . For 
example, without the knowledge of the fifth image, the 
innermost constraint given by the arcs on the density 
profile is the length of the radial arc that extends down 
to 3 arc-sec from the center. Its length depends on the 
source size which lies inside the caustic and needs to be 
related to the shape of its counter- image A5. A simple 
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Taylor expansion of the lens equation around the radial 
critical radius (where d rr <j> = 1) relates the half-length 
£ — 1.8" of the radial arc to corresponding source length 
du. This latter can be related to the size of the arclet A5 
(AA5 r = 0.5", AA5 t = 0.8") 14 which is triply imaged: 

du = ^-(1 — n r — r r n' r ) (B.7a) 
du = (1 — 2k c + R c + e c (3£ c + r c n' r — n c )) AA5 r 
+esm(2( V - <p ))(3£ c - 2n L ) AAb t . (B.7b) 



(|B.7a|l uses the property d rrr cj)(r r ) 



(1 — k t ■ 
V 7i)dx r 



r r K,' r )/r r whereas Eq. l|B.7b|) uses du = (1 — k- 
72 dxt ■ The equality between these two relations and the 

consti- 



normalization at the Einstein Radius (Eq. l|B.4a|P 
tute one more relation which is plotted on FigEl 

Appendix C: Adding haloes of galaxies as 
substructures 



In order to demonstrate that adding galactic halos to 
the models has weak impact on our conclusions regarding 
the cluster mass profile, we select the 9 brightest galax- 
ies which are the closest from the cluster center. Their 
I-band luminosities range between 0.17 < 10 -10 L//Lq < 
2.45 (the faintest has L w 0.17L*). We adopt a Faber- 
Jackson scaling to derive their respective velocity disper- 
sion. Each galaxy halo density profile is modeled by a trun- 
cated SIS with a cut-off radius r t . Without this trunca- 
tion, halos perturbations are constant and propagate up 
to the infinity. We adopt th e scaling laws proposed by 
ijNataraian et al.Smai i l2002t to study the lensing cluster 
AC114: 



K g( r ) = o 



6 = M- 



1/2 



1 



1 



\A- 2 +■ 



(C.la) 
(C.lb) 



with b related to a v as in Eq. I|13|l ■ ^* and rt* are two new 
free parameters. With this parameterization, perturbing 
galaxies have an individual mass-to-light ratio that does 
not depends on their luminosity. In the following, we only 
show the effect of perturbations on the softened isother- 
mal model, but we found similar conclusions for the NFW 
model. Fig. lC.ll shows the % 2 contours for this new couple 
of parameters after marginalization over the "macro" clus- 
ter model parameters whereas the constraints are the ones 
used in Sect. 13.4.21 and include the fifth image brightest 
peak knowledge. 

This modeling leads to a best fit x 2 /dof much closer to 
1. It shows also that introducing galaxy halo perturbations 
(in a way which is consistent with the radial, tangential 
arcs and their counter-images ) still predicts the fifth im- 
age at the observed position. The NFW case is similar. 

14 Note that the ellipticity of the mass distribution implies 
that the magnification matrix is not diagonal. Hence, the radial 
and tangential lengths correspond to the radial arc length. 



We illustrate on Fig. IC.lb the effect on the fifth image 
equivalent ellipse with the fiducial models referred as 1, 2 
and 3 on lC.lfa and compare it to the unperturbed softened 
isothermal model predictions. Galaxies haloes change the 
position and the shape of the fifth image only if they are 
so massive that they also damage significantly the exter- 
nal arcs image reconstruction. For instance, it yields to a 
bad % 2 /dof ~ 60, in the third model case. 
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Fig. C.l. Upper panel : xVdof contour plot for the couple 
(rt*,ff*). The "macro" model is the softened IS. Naturally, 
the \ 2 minimum (~ 1.1) is reduced as compared to the un- 
perturbed previous model. The three peculiar couples 1, 2, 3 
predict the fifth images plotted on the lower panel. The only 
noticeable effect of substructures on the fifth image occurs for 
models with a very high \ 2 model that provide a bad fit for 
the external arcs. 



